Bernoulli 14(1), 2008, 1-13 
DOI: 10.3150/07-BEJ6173 
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In this work, we give a closed form and a recurrence relation for a family of time-space harmonic 
polynomials relative to a Levy process. We also state the relationship with the Kailath-Segall 
(orthogonal) polynomials associated to the process. 
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1. Introduction 

Given a stochastic process X = {X t , t 6 R + } with finite moments of convenient order, 
a time-space harmonic polynomial relative to X is a polynomial Q(x,t) such that the 
process M t = Q(X t ,t) is a martingale with respect to the filtration associated with X. 
Major examples are the Hermite polynomials relative to a Brownian motion, the Charlicr 
polynomials relative to a Poisson process and the Laguerre polynomials relative to a 
Gamma process; for these and other examples, see Feinsilver [2], Schoutens and Teugels 
[11], Schoutens [10] and Barrieu and Schoutens [1]. 

Finding a family of time-space harmonic polynomials relative to a particular process is 
not an easy task; see, for example, Schoutens [10]. Nevertheless, when we are dealing with 
a Levy process X that has moment generating function in a neighborhood of the origin, 
Sengupta [13] mentions following Neveu [7] a general procedure based on the associated 
exponential martingale (the left-hand side of (1) below). Specifically, that martingale is 
an analytic function (of u) in some neighborhood of the origin and the corresponding 
Taylor expansion 

cxp{uX t } _ ^ u n 

in J y n ~ / , R n{Xt,t)— (1 

E[exp{uX t }\ ^ nl 

has coefficients such that R n (x,t) are time-space harmonic polynomials. However, when 
the process does not have a finite moment generating function, this approach does not 
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work; then, Sengupta [13] uses a discretization procedure to extend the results proven 
by Goswami and Sengupta [3], without explicitly showing the polynomials. 

In this paper, using Ito's formula in our proof, we give a closed form for a family of time- 
space harmonic polynomials of a Levy process with finite moments of all orders. From 
an interesting property observed by Sengupta [13], it is deduced that these polynomials 
play the role of the basic time-space harmonic polynomials because they span the set 
of all time-space harmonic polynomials. These results provide new insight into time- 
space harmonic polynomials and offer a general method for computing such polynomials. 
Furthermore, the construction has a direct extension to the case where the underlying 
Levy process has only a finite number of finite moments and we are only interested in 
a few space-time harmonic polynomials. We also present some examples and study the 
relationship between time-space harmonic polynomials and the orthogonal polynomials 
called Kailath-Segall polynomials (see Segall and Kailath [12] or Meyer [6]). 

The time-space harmonic polynomials that we propose here are related to the polyno- 
mials that give the moments of a random variable as functions of the cumulants. So wc 
will start with some results from the elementary theory of cumulants. 

2. Moments and cumulants 

Consider a random variable Y with finite moments of all orders. Denote by \x n (resp. K n ) 
its moment (resp. its cumulant) of order n. It is well known that fi n can be written as a 
polynomial in n\, . . . , n n , for example, 



If Y has moment generating function in some open interval containing 0, a general 
expression for fj, n can be deduced from the relationship between the moment generating 
function and the cumulant generating function: 



When Y does not fulfill the condition in the moment generating function, we can consider 
the finite series up to the convenient order; sec McCullagh [5], page 26, for the justification 
of this procedure. 

Through a direct examination of (2), Kendall and Stuart ([4], formula (3.33)) deduce 
that 



Mi = Kl, 



fl 2 = Ki + K2, 



fi 3 = n\ + 3kik 2 + k 3 , . . . 







where the second summation is over all non-negative integers r± , . . . , r m such that 

Em 
j=x jrj =n. 
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Let T n (x\, . . . , x n ), n > 1, be the polynomial denned by the expression (3). That is, we 
have 

/in = r ra («!,..., K„). (4) 

Also, write To = 1. 

A recurrence formula follows from (2) for T n ([15], Proposition 5.1.7): 

r„+i (xi , . . . , Xn+l) = ^ ( . J r 3 -(a;i, . . . , Xj) Xn+l-j ■ (5) 
3=0 

Further, 

g^- = ^JT n ^ j {x 1 ,...,x n - j ), j = 1, ... ,71. (6) 

From the Taylor expansion of r n (xi + y,x%, . . . ,x n ) at y = 0, we obtain the following 
expression that we will need later: 

r„(xi + y,x 2 , ■ . .,£„) = ^ f . \r n ..j(xi,. ■ .,x n _j)y J . (7) 
Interchanging the roles of x\ and y and evaluating the function at 0, we obtain 

r n (xi,x 2 , . . . ,x n ) = y~] ( . J r ri _j(o )x{. (8) 

3=0 

Applying this formula to Ki, . . . , K n and using the fact that k%, . . . , Kn are the cumulants 
of both Y and Y — EY~ gives the expression for the non-centered moment /t„ = E[F n ] (on 
the left) from the centered moments E[(Y — EF) ra_J ] and the powers of the expectation 
(EYy . It is a classical formula (see [4], formula (3.6)). 

3. Time— space harmonic polynomials for a Levy 
process 

Let X = {X t , t > 0} be a Levy process (meaning that X has stationary and independent 
increments and is continuous in probability and that Xq = 0), cadlag and centered, with 
moments of all orders. Note that Senguta [13] does not assume the stationarity of the 
increments in the definition of a Levy process. Denote by a 2 the variance of the Gaussian 
part of X and by v its Levy measure. For background on all these notions, we refer to 
Sato [9] . The existence of moments of all orders of X t implies that the Levy measure v 
has moments of all orders > 2. Write 

x n v(dx) for n > 2. 
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Observe that, for n > 3, m n is the cumulant of X\ of order n and m-i + a 2 is the cumulant 
of order 2. 

Following Nualart and Schoutens [8], we introduce the square- integrable martingales 
(and Levy processes) called Teugels martingales, related to the powers of the jumps of 
the process: 

Y t ^=X t , 

Y t [n] = (AX t )" - m n t, n>2, 

0<s<t 

where AX S = X s — X s -. The predictable quadratic variation of FW is 

(Y^) t = C n t, 

where C n is a constant. 

In the next theorem, we define a time-space harmonic polynomial related to the Levy 
process X. 

Theorem 1. Let X be a centered Levy process with finite moments of all orders. Then, 
for every n> 1, the process 

Mj = T n (X t , -(m-z + a 2 )t, -m 3 i, . . . , -m„i) 

is a martingale. 

Proof. Let n > 2. By Ito's formula, using (6) and the fact that [X, X]1 = a 2 t, we have 
M t (n) = nj* Mt~ l) dX s - Q (ma + a 2 ) J* ds 

- E (j) m i f M * n ~ j) ds + \ n ( n - 1 ) (j2 J* M * n ~ 2) ds 

+ ^2 (^n{X s - + AX S , -{mi + o- 2 )s,-m 3 s, . . . ,—m n s) (9) 

0<s<t 

- r n (X,_, -(ra 2 + cr 2 )s, -m 3 s, . . . , -m„s) 

- nAX s T n -i(X s -, -(ma + cr 2 )s, -m 3 s, . . . , -m„s)). 

Applying (7) and rearranging terms, we obtain 



T n (X s - + AI s ,-(m 2 + cr 2 )s, -m 3 s, . . . , - 




Time-space harmonic polynomials 5 
Then, the jumps part given in the expression for is 

E E ( n W- -i) (Ax s y = E ( n ) f^dxw 

0<s<tj=2 j=2 VJ/ JO 

=EQ^Mir j )(dy«+m idS ). 

J— 2 

Therefore, after some cancellations, (9) can be written as 

M (") = gQ^ M ir j) dF«. (10) 

Finally, since, for every t > 0, we < oo for k > and j > 1, all 

the stochastic integrals on the right-hand side of (10) are martingales. □ 

Corollary 2. Let X be a centered Levy process with finite moments of all orders. Define 
the polynomials 

Q n {x,t) = T n (x,-(m 2 + a 2 )t,-m 3 t, . . . ,-m„i), n> 1, 

and Qo(x,t) = 1. Then, 

(a) Q n (x,t), n> 1, are time-space harmonic polynomials relative to X; 

(b) Q n (x,t), n> 1, satisfy the recurrence relation 

Q n+ i(x,t) = xQ n (x,t) - n(m 2 + (T 2 )tQ n _i(x,t) 

~E ( •) m j+i*Q»-jOM); 

(c) an explicit expression for as a polynomial in X t is 
M t (n) =Q n (X t ,t) 



(11) 



E ( ^ ) r n-j(°> _ ( m 2 + cr 2 )i,-m 3 t, . . .,-m n _jt)X 3 t , n > 1; 



(d) z/ .X" is non-degenerate, then the family {Q n (x,t), n > 1} spans the set of all 
space-time harmonic polynomials relative to X . 

Proof. The recurrence in (b) is deduced from (5). Property (c) is obtained from 
(8). Finally, for (d), note that, from the properties stated in Section 2, the family 
{Q n (x,t), n > 1} satisfies the properties (i) to (iv) of Sengupta [13], page 953. On the 
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other hand, if X is non-degenerate, then the support of X t is unbounded (see [9], Theo- 
rem 24.3). Therefore the proof of Theorem 6 from Sengupta [13] can be applied to this 
family. □ 

Remark 3. 1. Assume that X is non-degenerate. Let X™ +p n -i(t)X"~ 1 + ■ ■ ■ +po(t) be 
a time-space harmonic polynomial relative to X, where Pj(t) are ordinary polynomials 
in t. Formula (2.2) of Goswami and Sengupta [3] can be transferred to our context. Then, 
the coefficients of the polynomial must satisfy 

Pi( s ) = X]( ■jPj( t )f J 'j-i( t ~ s )> 0<i<n, 0<s<t, 
j=i 

where /i r (i) = E[X t r ]. Applying this property to Q n (x,t), from (11) and the fact that 
fx r (t) = r r (0, (to2 + u 2 )t, m^t, . . . , m r t), it follows that 

r fc (0, -(m 2 + cr 2 )s, -m 3 s, . . . , -m k s) 

= ^2 ( \ ) Tfc-f (0, -(m 2 + a 2 )t, -m 3 t, . . . , —mu-it) 

x 1^(0, (m 2 + cr 2 )(t - s),m 3 (t — s), . . .,m e (t - s)), 0<s<t. 

2. When the Levy process X has only a finite number of finite moments, say up to 
order k, then Theorem 1 is still true up to the polynomial of degree [k/2] + 1, where [a] 
denotes the integer part of a. 



4. Examples 

4.1. Brownian motion 

Let {Wt, t>} be a standard Brownian motion and consider X t = Wt- Then, a = 1 and 
v = 0, so m n = for all n > 2. The time-space harmonic polynomials are 

Q„(x, t) = T n (x, -t, 0, . . . , 0), 

which satisfy the recurrence relation 

Q n+ i(x,t) = xQ n (x,t) - nt<2„_i(.T,i), n > 1. 

It is then clear that they coincide with the monic Hermite polynomials H n (x,t) (we 
follow the notation of Schoutens [10], Chapter 5) defined by 

cxp{ux — tu 2 12] = y H n (x,t) — - 
^— ' n\ 

n=0 



Time-space harmonic polynomials 
or, equivalent ly, 

H n (s,t)=^J H n {x/V2t), 
where H n (t) is the classical Hermite polynomial given by 

exp{2ux — u 2 } = \^ H n (x) — - 



n=0 



4.2. Poisson process 

Let {N t , t > 0} be a Poisson process of intensity 1 and X t — N t — t the compensated 
process. The Levy measure of X is a Dirac delta measure at point 1, hence m n = 1, 
n > 2. Then, 

Q n [X)tj (3*, £ , . . . , £) , 

and the recurrence relation is 

Q n+1 (x,t) = xQ n (x,t) -t^2 l n .)Qn-j{x,t). 

j= i \J '/ 

Now, consider the monic Charlier polynomials C n {x,t) (also with Schoutens' notation 
[10], Chapter 5), defined by 



exj n 

e-« f (l + «) I = VC B (x,i)^ 



It is well known that C n (N t ,t) = C n (X t +t,t) is a martingale. Write 

C n (x,t) = ^rC n (x + t,i). 

Then (see, e.g., [10], Theorem 8), 

C n {X s _,s)dX s = C n+1 {X t ,t). (12) 



Corollary 2(d) guarantees that the polynomial C n (x,t) is a linear combination of the 
polynomials Qj(x,t). To specify the relationship between both families of polynomials, 

we shall recursively construct the numbers A^™\ . . . , A„ such that 

n 

Q n (x,t)=J2^ n) Cj(x,t). 
3=1 
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First, 

Qi(x,i) =x = Ci(x,i), 

so A^ = 1. Assume that A^ , A^ , ■ ■ . , Aj™ 1 ,...,A^L 1 1 are known. Then, from 
(10), 



For the compensated Poisson process, all the Teugels martingales are equal, that is Y t 
X t , and using (12), we have 



(n) 



n—1 / \ n—j 

n 



< 1 , 

j = l ^ J ' k=l 



n—k 

' n 



k=i j=i 



It follows that A^ = 1 and 



n-1 

' n 



j=k 



4.3. Sum of two independent Levy processes 

Let {Y t , t > 0} and t > 0} be two independent centered Levy processes with mo- 
ments of all orders, and put X = Y + Z. Since the cumulant of order j of a sum of two 
independent random variables is the sum of the respectives cumulants of that order, the 
following convolution-type formula is deduced from (4): 

t n [Vl ' • • • ) Vn ~r Z n 

3=0 

Hence, 

where x = y + z. 

This result covers, for example, the process Xt = Wt + Nt — t, where {Wt, t > 0} is 
a standard Brownian motion and {N t , t > 0} is a Poisson process of intensity 1, both 
processes being independent. The polynomials corresponding to X are then convolutions 
between the Hermite polynomials and the polynomials studied in the previous example. 
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Let {Gt, t > 0} be a Gamma process, that is, a Levy process such that Gt has distribution 
Gamma with mean t and scale parameter equal to 1. Consider X t — Gt — t. The Levy 
measure is 

u(dx) = l( ,oo)(x)da; 
and the cumulants of X\ of order n>2 are 

/■OO 

m n = / x n - 1 e- x Ax={n-l)\. 
Jo 

The polynomials are then 

Q n (x,t) = r„(x, -t,...,-(n—l)\t) 
and the recurrence relation is 

n— 1 ^ 

Q n+ i(x,i) =xQ n (x,t) -n\t^2 — Qj(x,t). 

3=0 lh 

We are going to prove, in agreement with the results of Schoutens and Teugels [11] (see 
also [10]), that Q n (x,t) is related to the Laguerre polynomials. First, we compute the 
generating function of the polynomials Q n (x,t). Applying the general formula (2) to this 
specific case, for a in a neighborhood of zero of radius p € (0, 1), 

VQ„(i,i)^=Vr n (j;,-i,...,-( n -l)!^ 
* — ' n\ *■ — ' n! 

n=0 n=0 




= exp{u(x + 1) + 1 log(l — u)} 
= (1 — uf cxp{u(x + t)}. 

Second, using the relationship between Laguerre and Charlier polynomials, the generating 
function of the Laguerre polynomials can be written as (see [10] , page 47) 

oc 

(1 - uY cxp{ux} = Y^(-l) n L^- n \x)u n . 

n=0 

Then, 

Q n (x,t) = (-l) n n\L^- n \x + t). 
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Note that Schoutens [10] proves that L { t' n) {G t ) is a martingale, which is equivalent to 
the fact that Q n (X t ,t) is a martingale. 

4.5. Compound Poisson process with lognormal jumps 

The final example deals with a Levy process with finite moments of all orders, but 
without moment generating function. Let {N t , t > 0} be a Poisson process of intensity 
1 and {Y n , n> 1} a sequence of i.i.d. random variables with lognormal distribution of 
parameters fi = and a = 1, independent of the Poisson process. Consider the centered 
compound Poisson process 

3=1 

with the convention that the sum is when Nt = 0. As for all A > 0, E[cxp{A|Ai |}] = oo, 
it follows that X\ does not have finite moment generating function. The Levy measure 
is the law of Y\ and the cumulants of X\ of order > 2 are the moments of a lognormal 
law, 

m n = E[F"] = e™ 2 / 2 , n > 2. 

The polynomials are then 

Qni^^'i ^) — Ln(*£; ^6 3 • • • ) ^ ^ ) 

and the recurrence relation is 

Q n+1 (x, t) = xQ n (x, t)-tJ2 (™) c ^ +1)2 / 2 Q„_j (x,t). 

3=1 

5. Polynomials for a Levy process with moment 
generating function 

In this section, we study a general procedure (proposed by Sengupta [13]) for finding 
time-space harmonic polynomials related to a Levy process with finite moment generat- 
ing function in an open neighborhood of zero. We will prove that the resulting polynomials 
are the same as those obtained in Section 3. 

Assume that X\ has moment generating function in an open interval of the origin. 
Then, 

E[c uXt ] = cxp | iter 2 u 2 +t J (e ux ~ 1 - ua;)di/0) j 



Time-space harmonic polynomials 
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is analytic for u G (— <J, 5), for some <5 > 0, and it is known that (fixed u) 

exp{uX t } 
tW ' E[cxp{ M X t }] 

(13) 

= exp juX t - ^to- 2 u 2 - i J {c ux -1-ux) dv(x) 

is a martingale. On the other hand, for fixed u> € Q and t > 0, as a function of u, L t (u) 
is also analytic in (— <5, <5) and can be expanded as 



Lt(«) =53^(^,4) (14) 



Each coefficient R n (X t ,t) is a polynomial (in X t and t) that is a martingale (see the 
proof of a similar result in [10], page 46). Therefore, searching for a family of time- 
space harmonic polynomials is reduced to computing R n (x, t). This concludes the general 
procedure. However, we can go one step further since the condition that X\ has moment 
generating function in an open interval of the origin is equivalent (see Remark 4 below) 
to the condition that for some A > 0, 

exp{A|x|} v(dx) < oo. 

{H>i} 

Hence, the integral J R (e ux — 1 — ux) dis(x) is analytic in a neighborhood of and can be 
expanded as the power series X)^°=2 m n ■ Substituting the integral by the power series 
in (13), we obtain 

f 1 ^ u n 1 

L t (u) = exp I uX t - -tu 2 (o- 2 +m 2 ) - 1 2J ^j TO n \ ■ (15) 

Then, from (14), (15) and the exponential formula (2), we can identify the convenient 
/i„ and k n and obtain 

R n (X t ,t) = Y n {X tl -t(a 2 + m 2 ), -tm 3 , . . . , -tm n ) = Q n (X t ,t). 



Remark 4- The hypothesis that the random variable X\ has finite generating function 
in some open neighborhood of zero means that there exists A > such that 

E[cxp{A|Xi|}] <oo. (16) 

The function cxp{A|a;|} is submultiplicativc and therefore (16) is equivalent to 

/ exp{A|a;|} v{dx) < oo 

J{\x\>l} 
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(see [9], Theorem 25.3 and Proposition 25.4). It is worth noting that this last condition 
is equivalent to the hypothesis assumed by Nualart and Schoutens [8]. 



6. Time— space harmonic polynomials and orthogonal 
polynomials 

Another interesting family of polynomials relative to a Levy process X with moments of 
all orders is defined in the following way. Consider the iterated integrals 

P f (0) = 1, P t {1) =X t ,..., P t {n) = f Fir 1 ' dX s . (17) 

Jo 

These processes are related to the variations of X, 

xj»=X t , X^ = [X,X] t , X^=^2(AX s ) n , n>3, (18) 



s<t 



through the Kailath-Segall formula 



P t tn) = -{Pf-^XP - p^x[ 2) + ■■■ + (-l)™+ 1 P t (0) X t (n) ) (19) 

(see Scgall and Kailath [12] or Meyer [6]). 

It follows that pj n ^ is a polynomial in X^,...,X^ n ', is called the Kailath-Segall 
polynomial of order n. Denote this polynomial by P n (xi, ...,x n ), that is, 

P^=P n (X?\...,X^) 

and set Pq = 1. The translation of (19) to the polynomials P n (xi, . . . ,x n ) produces the 
recurrence relation 

Pn {x\ , . . . , Xji) 

(20) 

= -(P n -l(xi, . . . ,X n -i)xi - P n -2(XI,.. . 1 X n - 2 )x 2 H h (-l)" +1 Po^n)- 

n 

On the other hand, since the iterated integrals of different order of a Levy process are 
orthogonal, 

E[P t ^P t ^} = -C n t n 6 nm , (21) 

orthogonality of the polynomials P n (x[ 1 \ . . . , X^') follows. However, only when the 
underlying process is a Brownian motion or a Poisson process is P n (x[ , . . . , X t ) a 
polynomial of the form X™ + ai(t)X t n_1 + ■ ■ ■ + a n (t), with aj(t) being ordinary poly- 
nomials in t (see Sole and Utzet [14]). In light of (11), this is a major difference from 
time-space harmonic polynomials. 
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In general, both families of polynomials are related by the polynomial T n . Specifically, 

P n (xi,. ..,x n ) = —r n (x 1 ,-x 2 ,2lx 3 , . . ., (-l)" _1 (7i- l)\x n ) 
n\ 

because, by (5), this polynomial satisfies the recurrence relation (20). 
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